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We show how a fermionic quantum gas in an optical lattice and coupled to the field of an optical 
cavity can self-organize into a state in which the spontaneously emerging cavity field amplitude 
induces an artificial magnetic field. The fermions form either a chiral insulator or a chiral liquid 
carrying edge currents. The feedback mechanism via the cavity field enables robust and fast switch¬ 
ing of the edge currents and the cavity output can be employed for non-destructive measurements 
of the atomic dynamics. 
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The controlled generation of topologically non-trivial 
quantum phases is of greatest interest since they pos¬ 
sess special properties such as extended edge states that 
can be well protected against destructive environmental 
effects[1]. Thus, materials with topological non-trivial 
properties have the prospect to be utilized for quantum 
devices and quantum computing. Well known are topo¬ 
logical insulators with protected edge currents created in 
quantum Hall devices [2] by strong magnetic fields. 

For neutral atoms strong artificial magnetic fields can 
be created [3] which act similarly as magnetic fields for 
charged particles. Such artificial magnetic fields have 
been used in quantum gases confined to optical lattices 
to realize two showcase models of topologically insulating 
phases, the Hofstadter model in two-dimensions [4- 7] or 
on a ladder geometry [8] and the Haldane model [9]. 

Yet, the dynamic control and detection of topologically 
non-trivial quantum states remains a great challenge. In 
order to overcome this difficulty, in this work the dy¬ 
namical feedback between atoms and an optical cavity 
is employed to reach a self-organization of topologically 
non-trivial phases. One fascinating example for a self¬ 
organization of a coupled atom-cavity system has been 
realized recently by placing a bosonic quantum gas into 
an optical high-finesse resonator subjected to a perpen¬ 
dicular off-resonant pump beam [10-14]. Above a critical 
pump strength, the occupation of the cavity mode is sta¬ 
bilized and the bosonic atoms organize into a checker¬ 
board density pattern [12]. Many different proposals 
have been put forward to realize the self-organization 
of more complex quantum phases [13] reaching from the 
Mott-insulator [15] over fermionic phases [16-20] and dis¬ 
ordered structures [21-24] to phases with spin-orbit cou¬ 
pling [25-27], 

In this work, we engineer a direct coupling mechanism 
of the cavity photons to the tunneling of atoms in an op¬ 
tical lattice. This is achieved using a Raman transition 
induced by the combination of a transverse pump beam 
and the cavity field (Fig. 1). The frequencies of the pump 
beam and the cavity mode are chosen such that the en¬ 
ergy transfer is close to the energy offset between neigh¬ 
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FIG. 1. (color online) Sketch of the setup. Fermionic atoms in 
an optical cavity are subjected to an optical lattice potential 
which creates an array of ladders (see lower part) for which 
the tunneling amplitude along the legs is J«. The tunneling 
along the rungs is strongly suppressed by an potential offset 
A between neighboring wells. The tunneling along the rung of 
the ladders is restored by a Raman process using a transverse 
pump beam and the cavity mode (see upper part). 


boring lattice sites. The resulting process represents an 
effective tunneling of the atoms. Due to the running-wave 
nature of the pump beam, a phase f can be imprinted 
onto the tunneling of the atoms around a plaquette in 
the presence of cavity photons. We demonstrate that 
due to this novel coupling mechanism a self-organization 
into a state can take place in which the cavity mode is 
occupied and the atoms exhibit protected edge currents. 
Focusing on a situation in which the optical lattice is used 
to create an array of decoupled ladders (Fig. 1), either 
a chiral insulator or a chiral liquid are formed carrying 
chiral edge currents. Interestingly, the self-organization 
can set in depending on the filling either suddenly above 
a critical value of the pump intensity or can be exponen¬ 
tially activated with the pump intensity. In contrast to 
the application of two classical Raman beams [3], here 
the formation of the state relies on the dynamical feed- 
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back between the cavity field and the collective atomic 
tunneling. The cavity mode reaches rapidly its quasi¬ 
stationary state and by its extended nature could accel¬ 
erate the formation of correlations over the fermionic gas. 
Furthermore, the cavity output field can be employed for 
a non-destructive measurement of the tunneling of the 
atoms. 

Cold spin-polarized fermionic atoms are placed into 
a high-finesse optical cavity which is oriented along the 
^-direction (Fig. 1). We consider a single standing- 
wave cavity mode with frequency ui c , since the remain¬ 
ing modes are far detuned and do not contribute to the 
dynamics. Additionally, a strong lattice potential is ap¬ 
plied along the ^-direction which confines the atoms into 
decoupled planes. Within these planes, a superlattice 
structure is created using standing-wave laser beams with 
wavelength and 2X x along the x-clirection and X y along 
the y-direction. The phase difference of the two beams 
along the x-direction is chosen such that double wells 
with a potential offset A are created. The resulting ge¬ 
ometry is an array of decoupled ladder structures (Fig. 1) 
[8]. The atoms can tunnel along the leg direction of the 
ladder. In contrast, the tunneling process along the rungs 
is strongly suppressed by the potential offset between 
the two wells. To induce a cavity-assisted tunneling pro¬ 
cess along the rungs, an additional running-wave pump 
laser of wave length X p and frequency co p is applied par¬ 
allel to the legs of the ladders. This is detuned from 
the dispersively shifted cavity resonance frequency £j c by 
h(uj p — u> c ) ~ A and restores via a near-resonant Ra¬ 
man process the tunneling along the rungs of the ladder 
geometry (see Fig. 1). 

The cavity field and the pump laser are chosen suffi¬ 
ciently far-detuned from the atomic transition frequency 
w e such that the excited state is only weakly occupied 
and can be adiabatically eliminated. Thus, an effective 
description for the atoms in the internal ground state 
and the cavity mode can be derived (cf. Refs. [13, 28]). 
Performing additionally an expansion of the atomic field 
operator into Wannier states assuming that the optical 
lattices are chosen sufficiently strong and neglecting off- 
resonant terms, the resulting effective Hamiltonian (only 
a single ladder is considered) reads H = H c + H\\ + H ac 
with H c = h8 cp (Aa, 

H\\=—j\\ ( c in,j c ™,j+ 1 + h . c .^ 

0,1 

H ac = (aV^c^.c u + H.c.) 

j 

Here, c m j are the annihilation operators of an atom re¬ 
siding on leg m = 0,1 and site j. The term H\\ describes 
the tunneling of the atoms along the legs of the ladder 
with amplitude Jy. The operator a is the annihilation op¬ 
erator of a cavity photon in a frame rotating at frequency 
Lu p — A/h, and H c is the corresponding bare Hamiltonian. 


The two-photon detuning from the resonant tunneling 
transition between neighboring sites of a rung is denoted 
by S cp = [ w c — oj p + A /h] , and the AC-Stark shift caused 
by intra-cavity photons is neglected. The term H ac de¬ 
scribes the cavity-induced tunneling of the atoms along 
the rungs. Its amplitude is given by Htl = <j>j _, 

where f 1 p denotes the Rabi frequency of the pump and go 
the vacuum-Rabi frequency of the cavity. The factors <j>\\ 
and (f>± are effective parameters which can be tuned via 
the geometry of the optical lattice and the cavity mode 
[29], 

Due to the momentum transfer in the cavity-assisted 
Raman transition a site dependent phase <pj is trans¬ 
ferred to the atoms when they tunnel along the rungs. 
Thus, in the presence of a finite cavity field amplitude the 
fermions feel an artificial magnetic field oriented perpen¬ 
dicular to the x-y plane. The amplitude of this magnetic 
field is determined by the wave vector of the running- 
wave pump beam. For definiteness we choose X p = 2X y 
which results in a flux of p = 7r/2 enclosed by the plaque- 
ttes of the ladders. Other values of ip are also realizable. 
In contrast to previous setups [4-7], here the artificial 
magnetic field is not a static field imprinted by external 
laser beams, but dynamically emerges due to off-resonant 
scattering of pump light into the cavity mode. The re¬ 
sulting feedback mechanism between the rung tunneling 
process and the cavity field population leads to a self¬ 
organization of the system into a state in which the atoms 
are subjected to a strong artificial magnetic field. We will 
analyze the arising states in more detail in the following 
using a mean-field description for the cavity mode. 

Assuming that the cavity field reaches its steady state 
value rapidly compared to the timescale of the atomic 
motion, the expectation value a = (a) can be approxi¬ 
mated by the expression a = n _ i ^ (K±), where k is 

the damping rate of the cavity field. This means that on 
average the cavity field amplitude is proportional to the 
directed tunneling K± = e * v,?c o j c i j °f the fermions 
on all rungs. Note, that the complex phase of a can be 
chosen freely which reflects the 17(1)— symmetry of the 
problem. This will be spontaneously broken in an ex¬ 
periment, and in the following we choose without loss of 
generality (K±) > 0. 

The average value a of the cavity field can be sub¬ 
stituted into the equations of motion for the fermionic 
operators. The resulting fermionic evolution can be de¬ 
scribed by the effective Hamiltonian 

H f = H\\-(J ± K ± + K.c.). (1) 

The parameter J± = A(K\) with A = |P has to 
be determined self-consistently, because the expectation 
value of the directed tunneling (K±) does itself depend on 
the value of J±. A self-consistent non-zero solution only 
exists for a detuning S cp > 0, since only for a positive 
value of Ji the expectation value of the directed tun- 
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FIG. 2. (color online) Sketch of the band structure E± for 
different values of J± at ip = n /2. For J± < y2Ji\ the lower 
energy band has two minima Q± and depending on the filling 
two or four Fermi points occur. At J± = V2J\\ the energy 
bands separate. For J± > y/2 J \\, the lower energy band has 
one minimum at kd\\ =0 and a gap arises between the bands. 
The colored regions mark the Fermi sea for the filling n = 
1/2,1/4,1/8, respectively. Note, that these regions partially 
lie on top of each other. 


neling (K±) is positive. A non-vanishing solution of J± 
implies a finite expectation value of the photon operator 
by the relation a = J±(S cp + in)/(Ml5 cp ). The maximal 
feedback occurs for S cp = k because of the phase delay 
between the pump field and the scattered cavity field. 

For a fixed value of J± the effective fermionic 
Hamiltonian can be diagonalized [30-34], 

H F = T,k,a=± E <r( k )'Yl,k'Y<r,k with E ±( k ) 

—2Jy cos(</j/2) cos(M||)± \Jj± + 4J| sin 2 (</?/2) sin 2 (fcdy), 

where dy = X y /2 is the lattice spacing along the legs. 
The details of the Bogoliubov transformation are given 
in the supplemental material [29]. 

The energy dispersion for the case ip = tt /2 has two 
distinct regions (Fig. 2). For J± < \/2 Jy the two bands 
are overlapping and the lowest band has two minima at 

Q± = ±l/dy arccos (J + \ V For J± > \/2Jy 


the two bands are separated by a gap and the lowest 
band has one minimum at k = 0. 


One consequence is a transition between a liquid and a 
band insulator at half filling, n = N/(2L) = 1/2, where 
at J l = \[2 J|| (Fig. 2). Here N is the number of atoms 
and L the number of rungs on the ladder. At and be¬ 
low the filling n = 1/4 the occupied fermionic states lie 
entirely in the lowest energy band £L. 

The influence of the structure of the Fermi sea on the 
observables can be seen at the example of the rung tun¬ 
neling (Kj_) as shown in Fig. 3. At large J± V2J\\ 
and filling n < 1/2, the rung tunneling can be approx¬ 
imated by (K±)/L ss d||fcp/( 27 t), where ±fcp are the 
two Fermi points which depend linearly on the filling. 
At half filling the value of the rung tunneling increases 
slowly for small J± with a concave curvature until at 
J± = \/2Jy it shows a kink and crosses over to the lin¬ 
ear behavior. In contrast at n = 1/4 the rung tunnel- 


FIG. 3. (color online) The expectation value of the rung tun¬ 
neling (K±)/L versus the rung tunneling amplitude J± for 
different fillings n. The gray line represents the left-hand side 
of the self-consistency equation for a fixed value of AL. The 
crossings with the curves indicated by open symbols give the 
solutions to the self-consistency equation. 


ing shows a very smooth behavior in the entire range. 
It has a convex curvature at low J± given by a lin¬ 
ear rise with a logarithmic correction, i.e. (K±)/L « 


[log 32 — 2 log (Ji/J||)] 


Jl/J|| + o ( J±/J\\Y 


(2V 2 n) ---L/ -II T w L v--l/ -117 j ■ At lower fin¬ 

ing, e.g. n = 1/8, and low J± the behavior is similar to 
n = 1/4. However, at low values of J± the expectation 
value of the directed tunneling (K±)/L oc J±/J\\ and a 
slight kink occurs at the value Jj^o where the four Fermi 
points merge into two Fermi points. The value J± ^ de¬ 
pends on the filling and is given by J± o = dy for n = 1/8 
(Fig. 2). 

The gained insights into the dependence of the rung 
tunneling on J± enables us to find the solutions to 
the self-consistency equation. In particular, the self- 
consistency can be interpreted in the form J±/A = (K±), 
where the left-hand side stands for a linear curve with 
slope 1/A, whereas the right-hand side has the functional 
dependence of the rung tunneling. The crossings of these 
two curves give the solutions as indicated in Fig. 3. 

At half filling up to the critical pump strength of 
A cr L = 4-7T Jy/[v / 2/C(— 1) + /C(l/2)] R3 3.39Jy no non¬ 
zero solution of the self-consistency equation exists due 
to the concave curvature of the rung tunneling. Here, K, 
denotes the elliptic integral of the first kind. In contrast, 
for values larger than A cr two solutions exist over a cer¬ 
tain range (compare the crossings indicated in Fig. 3). 
The first solution corresponds to values J_l > and 

exists for all values of A > A cr . The second solution cor¬ 
responds to values J± < V2 Jy . This solution only exists 
in an intermediate regime of pump strengths. For both 
cases, the cavity field amplitude jumps to a finite value 
above the critical pump strength and the fermionic atoms 
are subjected to an artificial magnetic field (Fig. 4a). For 
the first solution the expectation value of the cavity field 
grows monotonically and depends approximately linearly 
on the pump strength, i.e. a ~ s n -i K ylliy l ar g e pump 
strength. 

The fermionic atoms form a band insulating state 
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FIG. 4. (color online) (a) The solutions J± of the self- 
consistency equation versus the pump strength AL for dif¬ 
ferent fillings n. Since the expectation value of the photon 
operator a is proportional to J_l, the curves represent also 
the dependence of the cavity field amplitude on the pump 
strength, (b) The induced chiral current J c versus the pump 
strength for different Filings. For n = 1/2 in both quantities 
a clear jump is present at a critical pump strength. The black 
and gray dashed lines show the two different solutions of the 
self-consistency equation in the case of half filling. 


which carries a chiral edge current which is de¬ 
fined as J c = Uj,o ~ where j jtm = 

—iJ\\ (cj] m Cj+i im — H.c.) is the local current on leg m. In 
contrast, for the second solution the expectation value of 
the cavity field decreases with increasing pump strength 
and the fermionic state corresponds to a liquid which 
carries a chiral edge current [35]. In both cases, a sud¬ 
den onset of the chiral edge current J c above the critical 
pump strength A cr occurs (Fig. 4b). At the jump the 
chiral current acquires its maximal value and decreases 
monotonically with increasing pump strength. 

In contrast, at quarter filling, n = 1/4, this abrupt self¬ 
organization changes into a smooth exponentially sup¬ 
pressed activation and the induced fermionic phase cor¬ 
responds to a chiral liquid. The self-consistency has for 
all pump strengths a non-zero solution which exhibits 
a slow initial rise, i.e. Jj_ = 4v / 2T||e -2v/ ^ 7r ‘ / ii/(" 4z d. At 
larger pump strengths it turns into the linear solution 
Jj_ ss AL d^k f / (2t:). Since the Fermi-momenta are lower 
than in the case of half filling, the reduced slope leads 
to a slower growth of the tunneling amplitude J±. As a 
consequence, the chiral current exhibits a slow rise with 
increasing pump strength (Fig. 4 b). 

At even lower filling, n = 1/8, only above 

a critical value of the pump strength A cr L ss 
y/2n Jy/log [cot(7r/16) tan(37r/16)] ~ 7.3Jy a non¬ 

vanishing solution of the self-consistency equation exists 
(Fig. 4 a). Above this critical value A cr a slow rise of J± 
occurs and a self-organization of the combined cavity- 
atom system into a chiral liquid takes place. The chi¬ 
ral edge current shows two different regimes of increase 
which are determined by the presence of two or four 


Fermi points and are connected by a slight kink (Fig. 4 
b). 

We briefly outline a possible experimental realization, 
using the fermionic isotope 40 K and cavity parameters 
k = 27rx4.5kHz and go = 2irx 0.76 MHz [36]. The lattice 
geometry is chosen such that Jy = 27r x 100 Hz, <j>» ss 1, 
(j>_ l > 0.05 and A/h = 2n x 23 kHz which is sufficiently 
large compared to the cavity line width. The pump wave¬ 
length A p = \ x is taken to be red-detuned by 1 nm from 
the D 2 -line at 767 nm and the cavity resonance such that 
S cp ~ 4re. Assuming an effective length of the ladder 
of L = 10 5 ( in the experiment several parallel ladder 
structures might be realized) the critical pump strength 
in the case of half filling, A CI L ss 3Jy, is reached for a 
pump Rabi frequency of 12“ = 27r x 100 MHz which cor- 

h 2 k 2 

responds to an AC Stark shift of 2AE r , where E r = 
and k x = 2 j n/\ x . Due to off-resonant electronic excita¬ 
tions, the pump light then causes per atom one spon¬ 
taneous scattering event every 0.7 s which is sufficiently 
large compared to the tunneling time h/J\\. Note that, 
since the effective coupling strength A scales with l/<5 2 p , 
the spontaneous scattering due to pump light can only 
be lowered by increasing either L, (f> j_, the cavity Purcell 
factor 7], or by decreasing Jy. 

Beyond the mean-field description, an effective dissi¬ 
pative dynamics with the jump operator K± with rate 
r ~ s 2 +k 2 cou ld drive the atomic system away from the 
ground state determined by Hp into a steady-state which 
is a dynamical equilibrium between driving and damping 
[13]. For the single-pump configuration considered above, 
this could (at long time-scales) lead - similar to optical 
pumping - to a transfer of the atomic population into 
the right leg of the ladders (see Fig. 1). This effect could 
be avoided by adding a second running-wave pump laser 
field along the y-direction which together with a second 
cavity mode (separated from the first cavity mode by 
twice the free spectral range) drives Raman transitions 
along the rungs into the opposite direction. 

The demonstrated self-organization of chiral insula¬ 
tors and chiral liquids in optical cavities could enable 
a dynamic control and non-destructive detection of such 
phases. In contrast to previous realizations, here, the ar¬ 
tificial magnetic field is generated by the feedback of the 
cavity mode. The quantum nature of the extended cav¬ 
ity field could lead to an effective (dissipative) long-range 
interaction between the fermions which could accelerate 
and stabilize the formation of the correlations. A gen¬ 
eralization of the discussed setup to a two-dimensional 
system [5, 6] and to interacting fermionic atoms can be 
envisaged and will be the subject of further investiga¬ 
tions. 

We acknowledge fruitful discussions with M. Fleis- 
chhauer, M. Kohl, H. Monien, H. Ritsch, and W. Zwerger 
and support from the DFG (INST 217/752-1 FUGG). 
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Bogolioubov transformation 

The Bogoliubov transformation is determined by the 
coefficients 



Coefficients from the tight-binding approximation 

In the limit of sufficiently strong optical lattice poten¬ 
tials, the field operators can be expanded into the cor¬ 
responding Wannier functions w y along the leg of the 
ladder and the maximally localized Wannier functions 
w m of the double well superlattice, where m = 0,1 labels 
the right or left well. To obtain the effective coupling 
strength which determines the amplitude of the tunnel¬ 
ing process along the cavity direction we can separate the 
contributions along the x- and (/-direction. Along the y- 
direction the following overlap integral is of importance 

h = J d ye~‘ lkvV w y {y)w* y {y). 

which is typically in the appropriate geometries close to 
one. Other contributions to the effective tunneling pro¬ 
cess along the x-direction from overlap intergrals along 
the y-direction can be much smaller and are therefore 
neglected. 

From the ^-direction two different contributions can 
play an important role depending on the exact geometry 
of the optical lattices [4-7]. The first one is given by the 
direct overlap integral between the Wannier functions of 
the two wells and the cavity mode function cos ( k c x ): 

= J dxcos (k c x) wi(x — d±)wo(x). (2) 

Here d± is the lattice spacing along the rungs of the lad¬ 
der. 

The second contribution is given by the time-periodic 
modulation of the energy offset between the sites of the 
double well [7]. Its amplitude is related to the on-site 
overlap integrals 

_ J J.j-gQg (k c (x + md±))w m (x)w’t n (x). 

Since in a double well potential with an energy offset, 
the values <fi and q are distinct, the coupling to the 
cavity mode induces a time-periodic modulation of the 
energy offset. Such a modulation leads to an effective 


tunneling along .x-clirection with its amplitude propor¬ 
tional to (/>^° — q to first order in the inverse of the 
oscillation frequency. 

In the main part of the article we use an effective pa¬ 
rameter 4 >_l which includes the two effects. The estimates 
given are minimal values which are determined from the 
overlap integral of the Wannier functions (Eq. 2). 
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